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We have extended our previous series studies of quantum antiferromagnets at zero temperature 
by computing the one-magnon dispersion curves and various structure factors for the hnear chain, 
square and simple cubic lattices. Many of these results are new; others are a substantial extension 
of previous work. These results are directly comparable with neutron scattering experiments and 
we make such comparisons where possible. 



is the archetypal model for describing long-range antiferromagnetic order in solids. Although there are no exact 
solutions in greater than one spatial dimension, a great deal is known about the model from various systematic 
approaches: exact diagonalizations, quantum Monte Carlo methods, and series expansions. Good overviews of the 
subject, with a particular focus on the square lattice and the relation to the high Tc cuprate superconductors, have 
been given by Barnes"'^ and Manousakis^. An area of particular current interest is the relation of models such as (1) to 
real materials. Quantities that can be most readily compared are the dispersion relations of low energy quasiparticle 
excitations and dynamical or integrated structure factors . The calculation of these is the main thrust of the current 
paper. At the same time the building of new and more powerful neutron scattering facilities is providing more precise 
data and allowing more detailed comparisons between experiment and theory'^''*'^. 

Our approach is through high-order 'linked cluster' series expansions^, where the quantities of interest are expanded 
perturbatively in powers of A (the so-called Ising expansion), and numerically evaluated at A = 1. This approach has 
been used with considerable success in computing ground state properties of quantum antiferromagnets^'*'^, and in 
computing the magnon excitation spectrum and spectral weight for the square lattice^*'. In our calculations we set 
J = 1 to determine the energy scale, except in comparison with experiment. In Section II we will define the various 
quantities of interest, and give a brief overview of the methodology. Section III gives new results for the structure 
factors for the linear chain. Section IV extends previous work for the square lattice*'^" and gives new results for the 
longitudinal and total structure factors. Section V gives results for the simple cubic lattice. Ground state series are 
extended by 2 terms and series results for the magnon energies and all structure factors are given for the first time. 
Finally in Section VI we summarize and attempt to relate our work to experiment. 



I. INTRODUCTION 



The spin- 



^ Heisenberg antiferromagnet, which we take in the exchange anisotropic form 




(1) 



II. METHODOLOGY AND DEFINITIONS 



The essence of the linked cluster method^ is the realization that many properties of a lattice model, in the thermo- 
dynamic limit Af — > oo, can be expressed as a sum of contributions from all possible connected or linked clusters of 
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sites which can be embedded in the particular lattice of interest. This is most obvious in the case of extensive bulk 
properties, such as the ground state energy, magnetization, susceptibility, etc, where we have 

Fr,{x)=Y,C{9/C)f,{x) (2) 
{g} 

where Fn{x) is the quantity of interest, with x representing the set of parameters in the Hamiltonian. The sum is over 
all clusters {g}, with C{g/C) being the embedding constant of cluster g in the lattice CoiN sites (proportional to N) 
and fg{x) a reduced quantity for cluster g. These latter quantities, which are independent of the lattice, are computed 
recursively^. It is easy to show that fg{x) is zero for any disconnected cluster, provided F is an extensive quantity. 
Linked cluster scries expansions arc then obtained by writing the Hamiltonian in the usual form for perturbation 
theory, H = Hq + AV^, and calculating the cluster contributions perturbatively, as series in A, up to some maximum 
achievable order (typically 10-20). The bulk series for Jjv(A) is then evaluated at fixed A, or extrapolated to A = 1, 
via standard numerical methods such as Pade approximants or integrated differential approximants^^. In practice all 
of this is done by computer and it is feasible to deal with of order 10^ distinct clusters. 

A stringent comparison between real materials and theoretical models is often provided by the spectrum of low 
energy excitations. These excitation energies can be measured in scattering experiments, and are characteristic of the 
quantum dynamics of the system. Gelfand^^ first showed how to compute excitation energies perturbatively, within a 
linked-cluster approach, and this is now a standard technique®. The basic idea is to compute an effective Hamiltonian 
matrix, which operates in the subspace of one-particle excitations of a cluster, use this to obtain a set of transition 
am,plitudes t{r) which describe propagation of the excitation through a distance r, obtain transition amplitudes for 
the bulk lattice by summing over clusters, and finally take the Fourier transform, giving the excitation energy in k 
space 

e(k) = X:i(r)e*''- (3) 

r 

While the dispersion relation (3) is an important probe of the quantum dynamics, an even more comprehensive probe 
is the dynamical structure factor 

1 f°° 

S^iK ^) = ^ y_ dte'^' E e*-'(5o (0)5?(i))o (4) 

i.e. the spatial and temporal Fourier transform of the dynamical spin-spin correlation function. The angular brackets 
denote an average (here a ground state expectation value), and a = x,y,z. This quantity is directly related to the 
cross section for inelastic neutron scattering (see e.g. Broholm & Aeppli^^). The integrated or static structure factor 
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'{SSS^) (5) 



is measured in an experiment where all neutron energies are included. 

For an isotropic system, in the absence of long-range magnetic order or other spontaneously broken symmetry, the 
components a — x,y, z oi (k, uj) or ^^(k) will be equal. This will no longer be the case if magnetic order is present. 
For a collinear ordered state, in the z direction, we need to distinguish between a longitudinal structure factor 

Si{k) = ^e*-[(5o^5^) - {S^){S^)] (6) 
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and a transverse structure factor 

St{k)^J2^"''"'(^SSr+sEsy) (7) 

r 

If unpolarized neutrons are used the cross section will measure the total structure factor 

Stotik) = Si{k) + StO^) (8) 

The dominant contribution to the transverse dynamical structure factor will come from one-magnon excitations, and 
St (k, Lo) will have the form 

St{k,w) = A,{k)6{uj - e(k)) + 5i„c(k,a;) (9) 

where j4i(k) is called the one-magnon spectral weight (or the exclusive structure factor) and S'incCk, w) is a smooth 
incoherent background term, arising from multi-magnon processes. It is easy to show that 

Ai(k) = l5^e*-(*o|(5o+ + 5o-)|*k)(*k|(5+ + 5-)|*o) (10) 

r 

where \^o), |^'k) are respectively the ground state and one-magnon state and S^, Sq are spin raising and lowering 
operators. It is also useful to define a relative multi-magnon spectral weight by 

m(k) = l-Ai(k)/St(k) (11) 

and a similar quantity for unpolarized neutron scattering 

Wtot(k) = l-Ai(k)/5tot(k) (12) 

The linked cluster formalism to compute the structure factor is relatively straightforward, and has been discussed 
in Refs. 14 and 10. The correlator sums 

Z^{r)^Y.(^rsr+r) (13) 

i 

are extensive quantities and thus have a linked-cluster expansion. There is, however, one interesting and important 
point regarding the longitudinal correlators and the structure factor. Linked cluster series for the correlators {S§S^), 
computed from a set of clusters up to some fixed maximum size, will have a maximum order in A which decreases 
with increasing r. On the other hand, the series for the compensated correlator {SqS^) — {Sq){S^) has a maximum 
order independent of r. This can be understood as follows. For any cluster the longitudinal correlator series all start 
with a constant (A°) term. Subtraction of subgraph contributions will cause cancellation of leading terms, leaving a 
series starting with some minimum power A^™'". However Pmin decreases with increasing r, and is zero for r = r„iax) 
the largest correlator which fits into the cluster, since, in this case, there are no subgraph subtractions. Thus, in the 
absence of the compensating term (14), much larger clusters would be required to give the largc-r correlator series to 
the same order. Inclusion of the compensating term avoids this problem since the leading terms in the bare correlator 
cancel and pmin (defined above), after subgraph subtraction, is independent of r. This allows longer series to be 
derived for the structure factor as defined in (6). The additional term 



^e^'^-(5o^)(5r^) (14) 
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will give a delta function peak at the antiferromagnetic wavevector Icaf, but will not change the longitudinal structure 
factor for k 7^ Icaf- The inclusion of this term reduces the total longitudinal structure factor, summed over momentum 
k, from to — A-P , where S and M are the spin and staggered magnetization, respectively. For the transverse 
structure factor this total sum is just S. 

There are two methods for computing series for the one-magnon spectral weight Ai(k). The first is to proceed 
directly from Eq. (10), as in Ref. 10. An alternative method^^ is from the linked cluster series for another quantity, 
the so-called 'exclusive matrix element', 

n{d) = (*o|(^+ + S-)\^m) ; S = ri-vm (15) 

where \^m) is the one-magnon wavefunction with initial unperturbed excitation at site m. Then 

Ai(k) = |^ri(5)e*-^|' (16) 
s 

The advantage of this second method is that it can be easily extended to the two-particle case, although we do not 
pursue this here. The two methods should, of course, result in the same final series. This provides a useful check 
on the correctness of the input cluster data, more stringent than the calculation of ground state bulk properties or 

excitation spectra. 

We compare our scries results with the prediction from spin-wave calculations. For the anisotropic Hamiltonian 
(1), the spin-wave theory has been computed to 4th order for the ground state energy, and 3rd order for most other 
properties'^. The second order spin- wave theory predicts the spin- wave excitation spectrum 

eic = zSq^, - I [C_igk + (A'^ - l)(C_i - C^){q-^ - q^)] (17) 

where z is the lattice coordination number, = (1 — A^7^)'/^, C„ is defined as 

Cn = ^Y.\^l->?llT"-l] (18) 



N ■ 

k 



and 



Ve-- (19) 



z 
p 

At A = 1, we can get a simple expression for the excitation spectrum 

ek = ^5(l-7^)'/'[l-C_i/(25)] (20) 

That is, the second order spin- wave theory only gives an overall a renormalization, with renormalization factor 
Zc = 1 — C-i/{2S), to the dispersion given by linear spin- wave theory. 
Linear spin-wave theory gives the transverse structure factor as 



In the limit fc = |k| — > 0, 7k — > 1 — k"^ /z, 



1-A 2A 



1/2 



so 5t(k) vanishes as Sk/^/2z at A = 1, while at k = 0, St vanishes as 5(1 — A)'/^/-\/2 as A ^ 1. 
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In the limit g = |q| = |kAF — k| — > 0, 7k — > — 1 + T 1^-, and 



1 - A , 2A 



1/2 



l/S.,k, = S-.(^^+,^^^J (23, 

SO S't(k) diverges as S^/2z/q at A = 1, while at k = kAF, 3% diverges as S^/2{l — A)~^/^ as A — > 1. 
We now turn to the series results. 



III. THE LINEAR CHAIN 

The anisotropic spin-^ Heisenberg antiferromagnet in one dimension (the XXZ chain) has been the subject of much 
study. Many materials which are well represented by this model have been identified (see Table 1 in Ref. 16). The 
possibility of exact results via Bethe ansatz methods has led to a good overall theoretical understanding of the model. 
In particular it is known that the elementary excitations are S = ^ spinons, or domain walls, with a dispersion 
relation^* 

espinon(fc) = /[cos2(fc) + 5^ sm''{k)]^^^ (24) 

where 

I = {l-Xy/'K{g'^)/n, 9'' = l-g' (25) 

and g is the solution of 

wK{g^)/K{g'^)= sech-\\) (26) 
and K denotes the complete elliptic integral, 

K{x)= / [l-xsm''{0)]-^/^de (27) 
Jo 

A series expansion for the spinon energy has already been derived by Singh^^, and shown to agree precisely with 
the expansion of the exact result (24) in powers of A. 

The structure factors are not known exactly for the XXZ chain, and here series expansions have a role to play. Singh 
et al}^ obtained long series for the longitudinal and transverse structure factors (6) and (7) at the antiferromagnetic 
wavevector A; = tt, to 22 and 12 terms respectively in A (only even terms occur in the longitudinal case) and studied 
the divergence of both quantities as A ^ 1 — . They found different exponents (^ 1.0,0.75) for the two power laws, 
and explained this apparently surprising result via a renormalization group argument. 

We have computed series for all of the structure factors, for general wavevector k, to order A^®. This represents 16 
additional terms for the transverse (and hence the total) structure factor. Our results for the isotropic case (A = 1) 
are shown in Figure 1. The structure factors diverge at fc = tt, as expected. For fc ^ tt, we find, to numerical accuracy, 
that Stot = 351, as expected, since the system has no long range order. For k = w, our longer series also show that 
longitudinal and transverse structure factors diverge with two different exponents, as found by Singh et al}^. 



IV. THE SQUARE LATTICE 

The square lattice S = ^ antiferromagnet has been much studied in recent years, largely due to its relevance to the 
high Tc cuprate superconductors. There is convincing, though not yet rigorous, evidence that the ground state has 
long-range Neel order, reduced by quantum fluctuations. 
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FIG. 1: The total and longitudinal structure factor for the linear chain. 
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FIG. 2: (Color online) The 1-magnon excitation spectrum e(k) along high-symmetry cuts through the Brillouin zone for the 
Heisenberg antiferromagnet on a square lattice. Also shown are the results of first order (blue dotted line), second order (red 
dashed line) and third order (green solid line) spin-wave theory. 

Some years ago we derived® perturbation series for the ground state energy, sublattice magnetization and parallel 
susceptibility to 14th order in the exchange anisotropic parameter A, and for the transverse (perpendicular) sus- 
ceptibility to order 13. These series provided very precise estimates of ground state properties for the entire range 
< A < 1, including the isotropic point A = 1. We also showed that higher order spin- wave theory was in excellent 
agreement with the series results. We have recently extended these series by two terms, to order A^^, the calculation 
involving a list of 185 690 clusters, up to 16 sites. We are happy to provide the new coefficients on request, but do 
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not present any new analysis of ground state properties here. 




FIG. 3: (Color online) The various integrated structure factors ^tot (unpolarized), St (transverse) and Si (longitudinal) along 
high-symmetry cuts through the Brillouin zone for the Heisenberg antiferromagnet on a square lattice. 

We have also extended an earlier calculation^" of the magnon excitation spectrum and spectral weight series by four 
terms, to order A^**. This calculation involves a large list of 4 654 284 clusters, up to 15 sites. The series coefficients 
are quite extensive and are not presented here, but we will provide them on request. We give in Table I the series 
at k = (7r,7r), (7r,0), and (7r/2,7r/2). The resulting magnon dispersion curve is shown in Fig. 2. It was obtained 
by extrapolating the series to A = 1, using integrated differential approximants. The first, second and third order 
spin-wave results®'^^ are included for comparison. We confirm the overall shape of the dispersion curve obtained 
previously^" but provide greater precision from the longer series. It is evident from the figure that the dispersion 
curve along the edge of the magnetic Brillouin zone (tt, 0) (7r/2, tt/2) is not fiat, as predicted by the first and second 
order spin-wave theory. We find numerically 

e(7r, 0) = 2.18(1) , e(7r/2, tt/2) = 2.385(1) (28) 

and so there is a 9.4% increase from (tt, 0) to (7r/2,7r/2). This agrees very well with a recent quantum A4onte 
Carlo calculation^" e(7i',0) — 2.16, e(7r/2,7r/2) = 2.39. Spin-wave theory, however, is unable to reproduce this 
variation even at third order^^ (via both Holstein-Primakoff and Dyson-Maleev transformations), which gives e{n, 0) — 
2.35858, e(7r/2,7r/2) = 2.39199. Our series results are also in qualitative agreement with experimental data for 
Cu(DC00)2 • 4D2O (CFTD)3 and Sr2Cu304Cl2^. However in La2Cu04 the observed magnon energy at (tt, 0) is 
higher than at (7r/2, 7r/2)^, opposite to the model result. It has been suggested^ that this is due to the presence of a 
significant ring exchange term in this material, but other explanations are possible^^. 

From our series for the magnon energies we can obtain a rather precise estimate of the spin wave velocity v. 
Following Singh and Gelfand^" we write the magnon energy at long wavelength in the form 



e{k) = C{X} + D{X)k^ + 0{k^) , k ^ \k\ ~* 



(29) 
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FIG. 4: (Color online) The 1-magnon spectral weight Ai, and multi-magnon spectral weights Wt and Wtot for the Heisenberg 
antiferromagnet on a square lattice. 

The spin wave velocity can be obtained from the series for 2C{\)D{X), evaluated at A = 1. Using integrated 
differential approximants^^, we estimate 2CD = 2.774(6) at A = 1, and conclude that v/Ja = 1.666(2). For k = 0, 
we expect the spin- wave energy to vanish as 



where the coefficient c can be estimated from our series: the result is c = 1.256(2). Third order spin-wave theory^^ 
gives v/Ja = 1.66802 and c = 1.23531, agreeing with the series estimates within 2%. 

In Figure 3 we show results for the various integrated structure factors along high symmetry lines in the Brillouin 
zone. The transverse structure factor was computed previously^" to order A^" - we have extended this series by four 
terms, to order A^^. Calculation of the longitudinal structure factor, and hence the total structure factor, by series 
methods is, as far as we know, given here for the first time. The series at k = (tt, tt), (tt, 0), and {tt/2, 7r/2) are listed 
in Table I. Various features deserve comment. Both longitudinal and transverse structure factors vanish at k = (0, 0). 
It is known, on general grounds, that the k dependence at this point is fc^, k respectively. Hence the longitudinal 
structure factor vanishes more rapidly. We estimate, from our series. 



where the coefficient for S't(k) is estimated using the same method as used for the spin- wave velocity v. A second-order 
spin-wave calculation^^ gives S't(k) = 0.10133fc. 

Both structure factors diverge at the antiferromagnetic wave vector k = {tt,tt). If the Neel state were an exact 
eigenstate the static longitudinal structure factor would be zero, except for a (5-function peak at (7r,7r). The actual 
shape reflects the additional contribution from quantum fluctuations. We first consider the asymptotic behaviour of 



e(k = 0) = c(l-A2)i/2, A^l 



(30) 



Si{k) = 0.042(4)fc2 as fc = |k| 
Stik) = 0.108(4)fc as fc |k| - 











(31) 
(32) 
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FIG. 5: (Color online) Ratio of longitudinal and transverse structure factors 251/51 for the Heisenberg antiferromagnet on a 
square lattice. Also shown, for comparison, is the QMC results^". 

longitudinal and transverse static structure factors at k = (tt, tt), as A ^ 1. Assuming 

5i(A) ~ (1 - A)-^' , 5t(A)^(l-A)-'^' (33) 

we estimate, from biased Dlog Pade approximants, that at — 0.50(2), while ai = 0.3(1). The exponents again differ, 
as in the ID case, but here it is ct; which is apparently smaller (this could be related to the fact that {S^) ^ on the 
square lattice). Linear spin- wave theory gives at — 1/2 (see Eq. (23)), but one would need a higher-order calculation 
to give ai, which has not yet been done. Next we consider the way in which the transverse and total structure factors 
at A = 1 diverge as k ^ (tt, tt). Defining q = (tt, tt) — k, we write 

St{q) = C{X) + D{X)q' + 0{ci^) , g=|q|^0 (34) 
Both C(A) and D{X) diverge at A = 1. However if we compute the inverse 

l/^t(q) = 1/C(A) - D{X)qyC^{X) + 0(q3) (35) 
and compare with the asymptotic form (see Eq. 23) 

l/^t(q)-[^(A)+i?(A)g2]i/2 (36) 

we find that S't diverges as {B^^^q)^^ with B — —2D/C^. The series for D for St is given in Table I. Our series, 
when analysed in this way, gives 

^t(g) =0.93(7)/g, q^O (37) 

The total structure factor series gives an estimate of 0.95(5), consistent with the same result. Spin- wave theory^^ 
gives 0.9288/?. 

Finally we note that the transverse structure factor exceeds the longitudinal one throughout the zone. The dominant 
one-magnon states only contribute to the transverse structure factor. The data can be analysed to extract the 1- 
magnon spectral weight Ai(k) and the relative multi-magnon spectral weights (Eqs. 11, 12). These are shown in 
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Figure 4, for the conventional lines in the Brillouin zone. The total 1-magnon spectral weight, summed over k, has 
the value 0.419(2), i.e. the 1-magnon excitations contribute 0.419/0.5 ~ 84% of the total transverse weight. We 
note that the maximum multi-magnon contribution to the structure factors, and hence to the integrated neutron 
scattering intensity, occurs at the (7r,0) point and is approximately 44% (29%) for unpolarized (polarized) neutrons. 
For k = (7r/2,7r/2), the multi magnon contribution is 31% (10%) for unpolarized (polarized) neutrons. Quantum 
Monte Carlo calculations'^^ give 40% (15%) at k = (tt, 0) (k = (7r/2, 7r/2)) for polarized neutrons. This is a significant 
contribution and needs to be allowed for in analysis of experimental data. 

In Figure 5 we plot the ratio 2Si/St throughout the zone. The overall shape is in excellent agreement with recent 
Quantum Monte Carlo data^°, but our maximum is about 0.62, considerably lower than the value 0.7 obtained by the 
Monte Carlo calculations^". Note that the Quantum Monte Carlo calculations have St/ Si diverging at k = (7r,7r), as 
they do not include the term (14) in their definition of the longitudinal structure factor. In principle, this term is a 
simple delta function at {tt,7t), and should not affect the measurement elsewhere for the bulk system. The omission 
of this term in the Monte Carlo calculations, however, can cause larger finite-size effects for finite systems, and this 
could be the cause of the discrepancy. 

V. THE SIMPLE CUBIC LATTICE 

We have carried out similar series calculations for the simple cubic lattice, and report on these here. Firstly, the 
previously calculated series for the ground state properties^ have been extended by two terms, to order A^**, involving 
a list of 180 252 clusters, up to 14 sites. This does not significantly change the previous estimates of ground state 
properties, and we do not present any further analysis. As usual, we are happy to provide the new coefficients to any 
interested reader. 
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FIG. 6: (Color online) The 1-magnon excitation spectrum e(k) along high-symmetry cuts through the Brillouin zone for the 
Heisenberg antiferromagnet on the simple cubic lattice. Also shown are the results of first order (blue dotted line) and second 
order (red dashed line) spin-wave theory. 
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Series for the magnon excitation spectrum have been derived, for the first time, to order A^". The calculations 
involve a list of 1487 597 clusters, up to 11 sites. The series for k — (7r,7r,7r), (tt, 0,0), (7r/2, 7r/2, 7r/2) are given in 
Table II. Figure 6 shows the magnon excitation spectrum along high-symmetry lines through the Brillouin zone, 
obtained from the series expansion, and first and second order spin-wave theory. It is evident from the figure that first 
order spin-wave theory gives the correct overall shape, but underestimates the magnitude by some 10%. The second 
order spin-wave theory is virtually indistinguishable from the series data, except on an enlarged scale along some cuts 
(as shown in the inset). A calculation of the spin- wave velocity, along the same lines as in the previous section, yields 
v/ Ja — 1.913(2). This compares with the first (second) order spin-wave value of 3^/^ — 1.732 (1.9003), and yields a 
quantum renormalization factor oi Z = 1.104(1) (compared to the square lattice with Z = 1.178(2)). This again is 
totally consistent with the lower relative effect of quantum fluctuations in higher dimensions. 
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FIG. 7: (Color online) The various integrated structure factors ^tot (unpolarized) , St (transverse) and Si (longitudinal) along 
high-symmetry cuts through the Brillouin zone for the Heisenberg antiferromagnet on a simple cubic lattice. 

Figure 7 gives our series estimates of the integrated structure factors along symmetry lines in the simple cubic Bril- 
louin zone. These are obtained from series expansions to order A^*'. We are unaware of any previous work along these 
lines. The same observations made for the square lattice can be made here. We note that the antiferromagnetic peak 
in Si is noticeably sharper here than for the square lattice, again reflecting the reduced role of quantum fluctuations. 
Finally in Figure 8 we show the 1-magnon spectral weight and the relative multi-magnon spectral weights. The latter 
are magnified by a factor of 10 for greater clarity. The multimagnon contribution to the transverse structure factor 
is nowhere greater than 3%, indicating the dominance of 1-magnon states, while the multimagnon contribution to 
the total structure factor, as would be measured by unpolarized neutrons, is as much as 15%. The total 1-magnon 
spectral weight, summed over k, has the value 0.482(1), i.e. the 1-magnon excitations contribute 96.4% of the total 
transverse weight. 

Similarly to the square lattice case, we obtain the following asymptotic results near k = and (tt, tt, tt): 



= 0.0114(2)fc^ k^ 



(38) 
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FIG. 8: (Color online) The 1-magnon spectral weight Ai, and multi-magnon spectral weights Wt and Wtot for the Heisenberg 
antiferromagnet on the simple cubic lattice. 




FIG. 9: Ratio of longitudinal and transverse structure factors 2Si/S't for the Heisenberg antiferromagnet on the simple cubic 
lattice. 

S't(k) = 0.1204(9)/;:, k^O (39) 
St{q) - 1.47(3)/g, q^O (40) 

Estimates from the 5tot series are consistent with these. 

In Figure 9 we plot the ratio 2Si/St throughout the zone. Here it has a maximum value about 0.3, substantially 
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FIG. 11: (Color online) Comparison of the 1-magnon transverse structure factor Ai(k) for CFTD^ (red solid points) and our 
series results. 
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FIG. 12: Comparison of the 1-magnon dispersion for Sr2Cu304Cl2^ and our series results with J = lO.SmeV, A — 0.976. 

VI. SUMMARY AND DISCUSSIONS 

The goal of this work has been to obtain numerically precise estimates of magnon energies and structure factors for 
the nearest-neighbour spin-i Heisenberg antiferromagnet for the linear chain {d = 1), square {d = 2) and simple cubic 
{d = 3) lattices. These quantities are directly comparable to experimental neutron scattering results, and the resulting 
comparison can provide a stringest test of the apphcabiUty of the simple model, as well as yielding an estimate of the 
(usually unknown) parameter J. 

We present such a comparison here for the quasi two-dimensional materials deuterated copper formate tetrahydrate 
(CuDC00)2 • 4D2O (CFTD)3 and the so-called "2342" compound Sr2Cu304Cl2.^ CFTD is a well characterized 2-d 
antiferromagnet"^. Figures 10 and 11 show a fit of our theoretical dispersion curve (Fig. 2) and 1-magnon transverse 
structure factor Ai(k) to the experimental data^, with a parameter J = 6.13meV. The overall agreement is very good, 
except near k = (tt, 0), where the theoretical one-magnon transverse structure factor is higher than the experimental 
results. The fitting parameter J is in good agreement with an earlier fit"^ to the previous series results"'^°. The strontium 
material is, a priori, more complex*. It contains two types of Cu^^ ions, Cui and Cun, and the interaction between 
these is fully frustrated. To the extent that one can regard these subsystems as decoupled, the Cun subsystem can 
be treated as an effective spin-i square lattice antiferromagnet with J ~ lOmeV. The measured dispersion curve 
shows a small spin gap, which can be modelled via a small magnetic anisotropy in the Hamiltonian. Figure 12 shows a 
comparison between the experimental data and our series results with J = lO.SmcV, A = 0.976, where A is determined 
from the minimum gap using Eq. (30). As is evident the fit is excellent, and again corroborates earlier results*. One 
should be cautious, however, in claiming too much from this and it would be highly desirable to have detailed structure 
factor data for further comparisons to be made. 

We are unaware of any good examples of spin-i antiferromagnetic materials with a simple cubic structure. 

Our results confirm, as expected, that the relative effect of quantum fluctuations decreases with increasing spa- 
tial dimension. Nevertheless, the multi-magnon contributions to integrated structure factors, and hence to neutron 
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scattering intensities, can still be appreciable even in three dimensions. 

For dimensions 2 and 3, the series expansion results arc in very good agreement with spin- wave theory, as far as 
it has been calculated. We conclude that the spin-wave calculations should be extended to higher order, to further 
check the agreement in quantities, such as the longitudinal structure factor, which have been little studied as yet. 

Note added: after this paper was submitted, we became aware of the work by Igarashi and Nagao^"^, who have 
performed a second-order spin-wave calculation of the transverse structure factor for the square lattice. 
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TABLE I: Series of square lattice one-magnon dispersion £(k), longitudinal structure factor ySi(k), transverse structure factor 
S't(k), and one-magnon exclusive structure factor Ai(k) at k = (tt, tt), (7r,0), (7r/2,7r/2), and series D for coeiRcient of (for 
e and S\ ) or (for St and Ai). Nonzero coefficients A" up to order n = 14 are listed. 



n 


k = (tt, tt) 


k=(77,0) 


k=(7r/2,7r/2) 


D 








dispersion e(k) 













2.000000000 




2.000000000 






2.000000000 




0.000000000 




2 


-1.666666667 




3.333333333x10"^ 






3.333333333x10" 


1 


1.000000000 




4 


3.171296296x10" 


1 


-9.953703704x10-2 






5.324074074x10" 


2 


2.569444444x10" 


1 


6 


-4.192337641x10" 


-1 


-1.693897891x10-3 






-9.073302469x10" 


-3 


6.581836259x10" 


1 


8 


2.709969904x10" 


1 


-2.806720342x10-2 






5.105325304x10" 


3 


—4.886280904x10" 


-2 


10 


-3.894335149x10" 


-1 


-1.062177000x10^2 






2.076320167x10" 


3 


7.984148319x10" 


1 


12 


4.289652578x10" 


1 


-9.046042545x10-3 






4.068538933x10" 


4 


-5.057247719x10" 


-1 


14 


-6.558882026x10" 


-1 


-8.872458632x10-3 






1.304340621x10" 


3 


1.460095528 










longitudinal structure factor S'i(k) 








2 


2.222222222x10" 


1 


1.111111111x10-^ 






1.111111111x10" 


1 


2.777777778x10" 


2 


4 


7.358024691x10" 


2 


-1.580246914x10-2 






3.703703704x10" 


4 


1.851851852x10" 


3 


6 


4.055166849x10" 


2 


1.167542552xl0"2 






-1.111176199x10" 


-3 


3.826627719x10" 


3 


8 


3.646524757x10" 


2 


1.065352379xl0"3 






4.750603835x10" 


3 


2.311003576x10" 


3 


10 


2.483688972x10" 


2 


2.298840789x10-3 






1.379013443x10" 


3 


1.874238376x10" 


3 


12 


2.211921434x10" 


2 


1.532346963x10-3 






1.422483566x10" 


3 


1.550055137x10" 


3 


14 


1.886279612x10" 


2 


1.141259073x10-3 






1.002657611x10" 


3 


1.314140042x10" 


3 








transverse structure factor S't(k) 











5 000000000x10" 


1 


5.000000000x10-^ 






5 nnnnnnnnn y 1 n" 


1 


000000000 




1 


6 666666667x10" 


1 


0.000000000 






0.000000000 




— 1 666666667x10" 


-1 


2 


3 333333333x10" 


1 


-1.111111111x10-^ 






-1.111111111x10" 


-1 


—2.222222222 x 10" 


-1 


3 


2.648148148x10" 


1 


0.000000000 






0.000000000 




—3 773148148 x 10" 


-1 


4 


2.390123457x10" 


1 


1.925925926x10-2 






1.086419753x10" 


2 


— 4. 340740741 X 10" 


-1 


5 


2.157488242x10" 


1 


0.000000000 






0.000000000 




—5.523841490 x 10" 


-1 


6 


1.922286407x10" 


1 


-1.297374724x10-2 






-6.735355253x10" 


-3 


—6.036591991x10" 


-1 


7 


1.825523316x10" 


1 


0.000000000 






0.000000000 




-7.022537598x10" 


-1 


8 


1.672647956x10" 


1 


-2.593884626x10-" 






-2.653241665x10" 


-3 


-7.488016779x10" 


-1 


9 


1.584816304x10" 


1 


0.000000000 






0.000000000 




-8.361027757x10" 


-1 


10 


1.488115597x10" 


1 


-2.819695062x10^3 






-1.579465770x10" 


-3 


-8.786542283x10" 


-1 


11 


1.436901069x10" 


1 


0.000000000 






0.000000000 




-9.576556162x10" 


-1 


12 


1.358134394x10" 


1 


-1.233261351x10-3 






-1.422867383x10" 


-3 


-9.970882006x10" 


-1 


13 


1.316662910x10" 


1 


0.000000000 






0.000000000 




-1.069785636 




14 


1.256571308x10" 


1 


-1.266236691x10-3 






-1.061011294x10" 


-3 


-1.106759542 










one-magnon spectral 


wcif 


;lit Ai{k) 











5 000000000x10" 


1 


5.000000000x10-^ 






5.000000000x10" 


1 


000000000 




1 


6 666666667x10" 


1 


0.000000000 






0.000000000 




— 1 666666667x10" 


-1 


2 


2 500000000x10" 


1 


-1.388888889x10-^ 






-1.388888889x10" 


-1 


— 1.944444444X 10" 


-1 


3 


1.425925926x10" 


1 


0.000000000 






0.000000000 




-3.245370370x10 


-1 


4 


3.326195988x10" 


1 


-1.143904321x10-2 






1.529706790x10" 


2 


-5.072723765x10" 


-1 


5 


3.917638154x10" 


1 


0.000000000 






0.000000000 




-6.858476264x10" 


-1 


6 


3.936459588x10" 


2 


-1.917201533x10^2 






-1.696167024x10" 


-2 


-4.528242195x10" 


-1 


7 


-8.018312217x10" 


-2 


0.000000000 






0.000000000 




-4.216442748x10" 


-1 


8 


3.778813233x10" 


1 


-1.855816875x10-2 






-4.208337928x10" 


-3 


-1.027026315 




9 


5.220328525x10" 


1 


0.000000000 






0.000000000 




-1.343102568 




10 


-1.799044627x10" 


-1 


-1.153727416x10^2 






-3.162792993x10" 


-3 


-3.701268973x10" 


-1 


11 


—d n.-^ysi.^/L/isvin- 


-1 


n nnnnnnnnn 






n nnnnnnnnn 




—d 71fi18Sn47vin" 


-2 
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TABLE II: Series of simple cubic lattice one-magnon dispersion e(k), longitudinal structure factor S'i(k), transverse structure 
factor ySt(k), and one-magnon exclusive structure factor j4i(k) at k = (tt, 7r,7r), (tt, 0, 0), (7r/2, 7r/2, 7r/2), and series D for 
coefficient of (for e and Si ) or (for St and A{). Nonzero coefficients A" up to order n = 10 are listed. 



n 


k = (7r,7r,7r) 


k=(77,0,0) 


k=(77/2,7r/2,7r/2) 


D 








dispersion e(k) 











3.000000000 




3.000000000 


3.000000000 




0.000000000 




2 


-1.950000000 




5.000000000x10"^ 


3.000000000x10" 


-1 


7.500000000x10" 


1 


4 


-7.480952381x10" 


-2 


5.483333333x10-2 


-1.550595238x10" 


-3 


2.859821429x10" 


1 


6 


-2.386949857x10" 


-1 


-1.194666672x10-2 


3.986473230x10" 


-3 


3.080615484x10" 


1 


8 


-3.884790029x10" 


-2 


2.654093113x10-^ 


6.579308964x10" 


-4 


1.842288801x10" 


1 


10 


-1.028725933x10" 


-1 


-5.327445001x10"" 


5.131262221x10" 


-4 


2.278014891x10" 


1 








longitudinal structure factor <S'i(k) 








2 


1.200000000x10" 


1 


4.000000000x10-2 


6.000000000x10" 


-2 


1.000000000x10" 


2 


4 


1.684444444x10" 


2 


-2.488888889x10-^ 


-3.962962963x10" 


-4 


1.530864198x10" 


4 


6 


1.301036907x10" 


2 


2.757538780x10-^ 


2.815062245x10" 


-3 


8.063018305x10" 


4 


8 


8.184237447x10" 


3 


7.116986056x10-"* 


1.126836657x10" 


-3 


3.866267816x10" 


4 


10 


6.005783346x10" 


3 


6.521285585x10"'* 


8.514828105x10" 


-4 


3.009882557x10" 


4 








transverse structure factor St (k) 











5 000000000x10" 


1 


5.000000000x10-^ 


5.000000000x10" 


-1 


000000000 




1 


6 000000000x10" 


1 


-2.000000000x10-^ 


0.000000000 




— 1 000000000x10" 


-1 


2 


3.000000000x10" 


1 


-2.000000000x10-^ 


-6.000000000x10" 


-2 


—1.200000000x10" 


-1 


3 


2.526666667x10" 


1 


4.614814815x10"^ 


0.000000000 




-1.887777778x10" 


-1 


4 


2.137481481x10" 


1 


-1.191111111x10-^ 


3.644444444x10" 


-3 


-2.033679012x10" 


-1 


5 


2.025150853x10" 


1 


-2.558718236x10-^ 


0.000000000 




-2.520957812x10" 


-1 


6 


1.752635491x10" 


1 


-7.923299511x10-" 


-3.718193643x10" 


-3 


-2.638276824x10" 


-1 


7 


1.685481230x10" 


1 


9.547987414x10-" 


0.000000000 




-3.039691386x10" 


-1 


8 


1.523090399x10" 


1 


-1.590182006x10-^ 


-9.298721272x10" 


-4 


-3.138723809x10" 


-1 


9 


1.480721342x10" 


1 


5.897338735x10-" 


0.000000000 




-3.487614896x10" 


-1 


10 


1.365201003x10" 


1 


-9.432640534x10-" 


-8.927606509x10" 


-4 


-3.574938809x10" 


-1 








one-magnon spectral weight At (k) 











5 000000000x10" 

KJ m\J\J\J\J\J\J\J\J\J y\ A.\J 


1 


5.000000000x10-^ 


5.000000000x10" 


-1 


000000000 




1 


6 000000000x10" 


1 


-2.000000000x10-* 


0.000000000 




— 1 OOOOOOOOOx 10" 


-1 


2 


2.812500000x10" 


1 


-2.875000000x10"^ 


-6.750000000x10" 


-2 


-1.162500000x10" 


-1 


3 


2.251666667x10" 


1 


5.309259259x10-2 


0.000000000 




-1.816944444x10" 


-1 


4 


2.268800324x10" 


1 


-1.446118552x10-2 


1.009873984x10" 


-3 


-2.105439590x10" 


-1 


5 


2.300857299x10" 


1 


-5.633309085x10-^ 


0.000000000 




-2.648295019x10" 


-1 


6 


1.606446880x10" 


1 


-2.164942522x10-^ 


-4.497506567x10" 


-3 


-2.550269987x10" 


-1 


7 


1.424582061x10" 


1 


2.355649892x10-^ 


0.000000000 




-2.873096589x10" 


-1 


8 


1.648269861x10" 


1 


-2.195243132x10-^ 


-1.493247410x10" 


-3 


-3.243370831x10" 


-1 


9 


1.707585298x10" 


1 


2.315207796x10-" 


0.000000000 




-3.678986537x10" 


-1 


10 


1.243609563x10" 


1 


-1.435173150x10-^ 


-1.241136381x10" 


-3 


-3.458041669x10" 


-1 




(7T,7r) (7T,0) (7t/2,7t/2) (tT.Tt) 



